Quantum Scattering in Quasi-ID Cylindrical Confinement 
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Finite size effects alter not only the energy levels of small systems, but can also lead to new 
effective interactions tmthin these systems. Here the problem of low energy quantum scattering by 
a spherically symmetric short range potential in the presence of a general cylindrical confinement 
is investigated. A Green's function formalism is developed which accounts for the full 3D nature 
of the scattering potential by incorporating all phase-shifts and their couplings. This quasi- ID 
geometry gives rise to scattering resonances and weakly localized states, whose binding energies 
and wavefunctions can be systematically calculated. Possible applications include e.g. impurity 
scattering in ballistic quasi-lD quantum wires in mesoscopic systems and in atomic matter wave 
guides. In the particular case of parabolic confinement, the present formalism can also be applied 
to pair collision processes such as two-body interactions. Weakly bound pairs and quasi-molecules 
induced by the confinement and having zero or higher orbital angular momentum can be predicted, 
such as p- and d-wave pairings. 

PACS numbers: 03.65.Nk, 05.30.Fk, 05.30.Jp, 34.10.-l-x 



I. INTRODUCTION 

Matter at very small dimensions can be strongly af- 
fected by boundary and surface effects. Bulk 3D prop- 
erties may then change substantially or even disappear. 
Reaching this regime poses a long term challenge moti- 
vated among others by the continuous technological drive 
towards ever smaller information processing devices 0. 
Below a certain limit of small length scales, novel quan- 
tum electronic properties may appear and be useful in 
order to develop alternative devices ■ 

Recent developments in atomic physics have addressed 
similar questions by using atom-optical devices (or atom- 
chips) made out of surface-patterned substrates 0, 0, 
IE IE 111 or by means of laser manipulation of atoms 
for trapping and atomic lithography ^1^. Analogous 
to electronic matter waves, these devices and techniques 
have the potential to trap and manipulate atomic matter 
waves down to sub-micron scales. A vast range of phe- 
nomena with atoms can then be studied, several of them 
bearing close relationships to specific solid state systems. 

Of particular interest for both bulk and atomic sys- 
tems is the quasi- ID regime of quantum coherent sin- 
gle mode transport. In this ultimate limit, the degrees 
of freedom in the transversal trapping directions are ef- 
fectively frozen in the quantum-limited ground state of 
the confining potential. Such one-dimensional system 
can show distinct physical properties e.g in quasi-lD 
quantum wires in two-dimensional electron gas (2DEG) 
systems H El E3JllJll_|ll EE E3 and Tomonaga- 
Luttinger liquids |l8l Il9l l2ol|~^In elongated gases of ul- 
tracold atoms, the one dimensionality is revealed e.g. in 
magnetic guides |EfllEJEJilE| , in phase fluctuations of 
quasi-condensates |2ll 123. |2J . t he Tonks-Girardeau gas 
of impenetrable bosons |2 4 l25l I M Wk I2I. |29| or one 
dimensional optical lattices |30l l3l| . In this broad con- 



text, we address two processes that are affected by the 
quasi- ID geometry. On one hand, the matter waves can 
be scattered off e.g. impurities or defects fixed along 
the guide. This scattering by central fields is an im- 
portant variable in the transport properties and can be 
strongly affected by the confinement. This has been re- 
vealed by several studies of mesoscopic systems in the 
context of the Landauer approach to the quantized con- 
ductance of quantum point-contacts in the presence of 
impurities 0, 0, 0, 0| . On the other hand, two-body 
collision processes should be dealt with when the mat- 
ter wave consists of interacting particles. A full under- 
standing of these low dimensional systems should there- 
fore deal with both types of scattering. It should also 
derive ad hoc phenomenological parameters such as effec- 
tive low dimensional coupling strengths from the real 3D 
character of the interactions. Another context in which 
such dimensional reduction plays a role is the transmis- 
sion and reflection of electromagnetic 32] and elastic 3^ 
waves in transmission lines or resonators. 

New light has recently been shed on this subject by 
studies of atomic collisions in ID |3 IsE HE Hzl and 
2D |2E 1^ JQ] ■ Focusing on a parabolic cylindrical con- 
finement, for which the center of mass can be eliminated, 
pioneering predictions for collision processes are possible, 
e.g. that of confinement induced scattering resonances 
(CIR) in ID. In such a resonance HE lUli a total 
reflection may take place between the colliding atoms. 
The s-wave (or zero-range) approximation to the scat- 
tering potential 's^, 'sE'l may be partially lifted, provided 
the coupling of orbital angular momenta due to the cylin- 
drical geometry are not considered [sEl- Another prop- 
erty induced by the conflnement and closely related to 
CIR is the formation of weakly bound quasi-molecules 
with zero (s-wave) orbital angular momentum |37| , which 
have been recently observed with singlet fermionic atoms 
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FIG. 1: (a) Longitudinal profile of the cylindrical guide. 
The scattering potential V{r) is limited to within the spher- 
ical shaded area at the center, (b) Total potential energy 
along the transverse a;y-plane passing through the origin r = 
(thick curve). For an attractive V{r), one deep Eo and an- 
other shallow Eb bound state without confinement are indi- 
cated, with their respective probability density profiles (dot- 
ted curves). 



in tight laser traps 31]. Coupling to the center of mass 
under non-parabohc confinement in the s-wave approx- 
imation is treated in 41] and can lead to the appear- 
ance of additional CIRs. In the independent and dis- 
tinct context of impurity scattering studies, rectangular 
geometries predominate which are most suitable for the 
growth of semiconductor heterostructures. Similar prop- 
erties of resonant scattering and weakly localized states 
are then predicted for non-interacting 2DEG systems. 
This is demonstrated in the s-wave approximation to the 
scattering potential in . Finite range scatterers under 
general rectangular confinement are also considered, but 
in an already reduced two-dimensional space [T^ Hsl IT^ . 

In this work, a detailed formalism is presented that 
extends the results for impurity scattering and collisions. 
It can treat general cylindrical confinement, parabolic or 
not, and incorporate all scattering phase-shifts beyond 
s-waves, as well as the full couplings of orbital angu- 
lar momenta due to the broken spherical symmetry. A 
comprehensive assessment of the scattering process re- 
veals the most important mechanisms and parameters at 
play. From this unified description, CIRs are seen to be 
a general low energy effect in quasi- ID geometries. It 
might be useful e.g. as an alternative gating mechanism 
in low power transistor-like devices that could incorpo- 
rate quasi- ID structures |l[ll[ll[ll[ll[ll[ll[Il|43- 
The existence of weakly localized bound states and quasi- 
molecules can, in turn, be extended to include higher or- 
bital angular momentum pairing, whose binding energies 
and wavefunctions can be systematically computed. 

This paper is organized as follows. In section ^1 the 
present problem of scattering in confined geometries is 
formulated for both scattering by central fields and colli- 
sion processes. In section Hill the s-wave approximation 
to the scattering potential is briefly reviewed. In sec- 
tion llVI after an overview of the free space scattering, our 
formalism in terms of Green's functions and phase-shifts 
is developed. In section this formalism is applied to 
analyze two types of confinement induced effects, namely. 



CIR and weakly localized states. A discussion of the main 
approximations used is given in section IVll followed by 
the conclusions. 



II. SCATTERING IN LOW DIMENSIONALITY 

In a low energy scattering process under confinement, 
one expects the scattering potential to affect primarily 
the degrees of freedom not experiencing external confine- 
ment. The other degrees of freedom are then forced to 
occupy only one or very few confined states. However, 
the very existence of these confined states do play a role 
due to virtual transitions. The effects of the latter can be- 
come noticeable when properly examining the unconfined 
degrees of freedom. In this regard one is lead to calculate 
effective scattering amplitudes as given by Ea. f5b|) . 

These amplitudes provide a complete description of the 
quasi-lD problem. For real energies, they describe the 
scattering properties of running waves. Weakly localized 
states are then found in the complex-energy plane. 

A. Hamiltonian under confinement 

Consider first the scattering of a single particle of mass 
/i and coordinates r = (x, y, z) by a central field fixed at 
the origin. Let V{r) be the spherically symmetric scatter- 
ing potential and U (p) the cylindrically symmetric con- 
fining potential, where p = {x, y) denotes the transverse 
coordinates. The Hamiltonian is then given by 

iJ = -^V2 + C/(p) + y(r). (1) 

If Rv and Rjj are the ranges of V{r) and U{p), respec- 
tively, then \V\ < Vb for r > Ry and t/ < C/q for 
p <C Ru , where Vq and are characteristic energy 
scales, such as the respective ground state energies (see 
fig. Two particular examples are the parabolic con- 
finement U{p) = p,uj-^p- /2, for which Ru is of the order 
of the oscillator length a± = {h/ pLu±y^^ and the square- 
well (or hard- wall) such that U{p) — > -l-oo for p > Rjj 
and zero otherwise. The first condition on the scattering 
and on the confining potentials assumed in this work is 
Rv ^ Ru ■ In this sense, both the potentials V{r) and 
U (p) can be relatively general for this type of scattering 
process. This condition implies that there is a distinct 
region Ry ^ r <C Ru where spherical symmetry still 
prevails and one can define scattering phase-shifts. This 
is an (intermediate) asymptotic region, where the effect 
of the confinement U{p) is felt only as a boundary condi- 
tion to the scattering by V{r) (see section llV C|l . Fig.^ 
depicts the geometry of the problem. Another important 
condition for the validity of the present approach to scat- 
tering in confined geometries is the condition k ~ 1/Ru 
of low total scattering energy E = h'^k'^/2p (or Fermi 
energy) for the Hamiltonian H. Only few excited trans- 
verse states (fin above the ground state (see Eq.(|2Jl) can be 
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effectively populated. The scattering process is thus as- 
sumed to occur at low temperatures or under sufficiently 
tight confinement. 

Alternatively, consider the case of collisions between 
e.g. two cold atoms of masses m and coordinates = 
{pi.Zi), z = 1, 2 . R = {ri + r2)/2 and r = ri - r2 
denote the center of mass (CM) and relative coordinates, 
respectively. The total Hamiltonian is 



V{r), 



where Uc is the cylindrical confinement and V the two- 
body interaction potential. For this scattering process 
Uc should be parabolic, so that the CM can be separated 
exactly. The Hamiltonian for the relative coordinates is 
then given by Eq.JQ), where fi = 171/2 becomes the (col- 
lisional) effective mass and U{p) — 2Uc{p/2) . For non- 
parabolic confinement, additional scattering resonances 
that can occur even at weak coupling to the CM is 
not included in Eq.(^. However, at ultracold tempera- 
tures for which only the lowest transversal states of the 
guide are occupied, the probability to find the atoms in 
the center of the guide is largest. Consequently, it might 
be a reasonable zeroth-order approximation to assume 
on average that the CM coordinates are also confined to 
the center of the guide. Therefore Eq.0J is a starting 
point to the dynamics of the relative coordinates for a 
general interaction V{r), in which the CM fluctuations 
around the mean value (R) = are at the present stage 
neglected. 

Because of the symmetry of U{p) + V{r), the scat- 
tering solutions can be restricted to be axially symmet- 
ric. The corresponding axially symmetric and orthonor- 
malized eigenstates of the confined variables ipn{p), n = 
0, 1, 2, . . . , satisfy 



ipn{p) ^ in iPn{p) (2) 



where e„ = h^q^/2ii > are the transversal energies, 
with qo ^ 1/ Rij . 



B. ID scattering 

In a purely ID problem let the ID interaction potential 
Vid{z) have the finite range Rid around the origin z — Q, 
such that ViD ~ for \z\ ^ Rid- The scattering problem 
for an incoming plane wave e''^°^ can be described by 
the amplitudes in the asymptotic behaviour of the 
solution 'ip{z) in the transmission and reflection regions 



i^{z) 



ikoz _|_ y-g-ifcoz 



Rid < z, 
z < -Rid- 



Current conservation then yields 

|i + /+P = i-irr 



(3) 



Particular cases are resonant transmission (/+ 0) and 
total reflection (/+ ^ — 1) at finite feg and finite poten- 
tials. 



C. Effective ID scattering 

The above ID picture corresponds to a fixed unper- 
turbed transversal state, e.g., the ground state ipo{p)- 
The total wave function would then have the form 
'ip{z)ifQ{p) where scattering effects would occur only in 
Tpiz)- However, transitions to other states (^„(p), n = 
1,2,..., even as virtual transitions, arise from the cou- 
pling between the 's caused by the scattering potential 
V{r). 

Let 5'(r) be the full 3D axially symmetric scattering 
solution to Eq.lQJ, 



[V^ - u{p) + fc2] «'(r) = v{r) *(r). 



(4a) 



where E = h^k^ /2p > is the total scattering energy, 
u{p) = 2/LtC/(p)/ft2 and v{r) = 2pV{r)/h'^. The above 
mentioned couplings and virtual transitions are best seen 
by expanding the solution in the cylindrical basis defined 
by the (/3„'s, 



X! i^n{z)iPn{p), 
n=0 



(4b) 



and solving for each ipniz). Substituting into Eq.(|4a| 
multiplying by (p*, integrating and using Eq.Q yields 



dz^ 



dx dyipl{p)v{r)'^{r) 



For a given total energy E, one defines the longitudinal 
wave vectors fc„ by 



if < n < riB , 
if n > ue-, 



where is the largest integer such that < k < 
qi+UE and denotes the highest "open channel" . In the 
single mode regime ue = 0- As usual, the general so- 
lution ipniz) can be expressed in terms of a homoge- 
neous solution and a particular inhomogeneous solution, 
Tpf^niz). This inhomogeneous solution can be written in 
terms of the ID Green's function 



G„(z) 



2'ikf^ 



+ k'^Gn — —S{z — z'), 



describing an outward scattering for n < and an expo- 
nentially decaying virtual state for the "closed channels" 
n> He- Then for all n 

i^Uz) = - 1 d'r' G^{z - z') vl{p') v{r') ^{r') . 
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FIG. 2: A pictorial view of a scattering "path" of a par- 
ticle coming from the left. Outside the effective ID range 
\z\ 3> R'idj the particle is stabilized by the confinement back 
into e.g. the tansversal ground state (for the single mode 
regime ue ~ 0). Here the net effect of the scattering is con- 
tained in the amplitudes and fg of backward and forward 
scattering, respectively. Within the dashed circle, however, 
the spherical symmetry prevails. The phase information of 
the wave function within this region is then captured by ff^. 



The expansion Eg . (|4b|l of ^'(r) takes then the form 



n=0 



^-'fAz)Vn{p) , (4c) 

n=0 



where the homogeneous part &„e*'^"^ of each ipn{z), for 
some constant 6„, gives rise to the total incoming state 
and thus is limited to the open channels n < tie only. 

The asymptotic condition is obtained in the limit l^l — > 
oo by neglecting the exponentially decaying terms ipf^niz) 
in Ea. (|^ for n > tie- Using then \z — z'\ — \z\^ z' in 
Gn, for z —^ ±oo, one has 



*(r) 



71 E 



n=0 



AknZ 



(5a) 



where the effective ID scattering amplitudes are de- 
fined by (for n <nE) 



the scattering center, such that Ea. (|5a() is the dominant 
contribution to the scattering wave function (see fig. El- 
As regards collision processes, this effective range sets 
some limits on the average equilibrium distance between 
the colliding particles. 

The present formalism to quantum scattering under 
confinement is approximate but non-perturbative, in the 
sense that there is no small parameter around which 
an expansion is performed. The determination of its 
range of validity is chosen here by imposing the conser- 
vation of probability on the final wavefunctions. The 
probability conservation condition in 3D and the iden- 
tity d'^r divj{r) = cPS ■ j{r) for a volume 
il enclosed by a large cylinder of surface area Sc and by 
two transversal discs of surface areas Si at zi < and 
52 at Z2 > yields 







/ d(t)dzp[p-j{r) 



dxdy k-j{r) + dxdy k-j{r) 

Si L J "'52 



where p and k are cylindrical unit vectors. Applying 
^I' = V'n'/'n, thc first integral on the rhs vanishes for 
large radius p — > oo since each (pnif) — > 0. The other 
two integrals yield 



= 



El 

ri=0 



'3n{zi) + jn{z2)] 



where jn{z) are ID currents calculated with ^n{z). Set- 
ting zi^2 ^ TR'id oii^ the asymptotics for -0„ in 
Ea. (|5a| so that the total conservation condition reads 



0^22{\bn + f+\'' + \f-\'-\bn\')k„ 



(5c) 



J n 



2z kfi 



d^i 



^„(p') i;(r')*(r'), (5b) 



for forward z > and backward z < scattering, re- 
spectively. Note first that although \z\ 3> \z'\, the phase 
knz' in is not necessarily negligible since kn\z'\ ^ 2tt 
for sufficiently large momentum kn and depending on the 
range Ry of the interaction v{r'). This phase is negligible 
only when the wavelength 27r/fc„ of the incoming wave is 
not able to resolve the details of the potential u(r') at low 
momenta. In this case, a zero-range or s-wave approx- 
imation to v{r') is valid. However, this approximation 
should break down when the scattering occurs at short 
wavelengths or when higher angular momenta become 
necessary. Note that there is no clearly defined ID range 
Rid as in the purely ID case, since the asymptotics given 
by Eq. H5a|l depends on the convergence behaviour of thc 
series in Ea. (|^ . If the convergence is slow, one expects 
an effective ID range, say R[ei substantially larger than 
Rv or even Ru. It is an effective range for the ID colli- 
sion and is of the order of the minimum separation from 



This equation extends the pure ID result in Eq.Q. It 
also determines how the initial flux is distributed among 
the open channels. In other words, even when the initial 
state occupies only a single channel, i.e., 6„ = So^n, new 
channels are occupied due to scattering, if the total en- 
ergy is large enough ue > 0. This conservation condition 
can serve then as a gauge for the range of validity of our 
approximations to the amplitudes f^. 



III. SCATTERING IN THE s-WAVE 
APPROXIMATION 

In many circumstances, the low energy scattering in 
free space can be well described by approximating the 
interaction V{r) by Fermi-Huang's zero-range poten- 
tial illii 



Vs{r) = <5(r)-^(r • ) 

jj, or 



(6a) 
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This approximation singles out the I = angular momen- 
tum component of the collision and thus depends only 
on the s-wave scattering length a. The regularization 
d{r ■ )/dr assures that the singularity 1/r of the full 
scattering solution to this Vs{r) is properly dealt with. 

For collisions in a strong parabolic confinement and 
using Eq. Ht)a|) , the series in Eq. also has a singularity 
for short distances r — > very similar to 1/r. By properly 
dealing with this singularity, with |34| or without any 
regularization [35|, the amplitudes in Ea. (|5b|l can be 
calculated. In the single mode regime {ue — Q) and for 
small longitudinal momenta fco one obtains 



1 



l + i 



2a ' ^ 



fco 



(6b) 



where C — 1.4603... and a±_ is the harmonic oscilla- 
tor length of the parabolic confinement [s^. l35l | . Here, 
a total reflection = —1 is predicted for large scat- 
tering lengths, on the order of the confinement length, 
such that aj_ — Ca , which is called the confinement in- 
duced resonance (CIR) 34]. This total reflection means 
that the colliding pair experiences a diverging effective 
ID interaction along the longitudinal cylindrical axis (see 
section IVD 

It will be seen later that such a singularitiy should ac- 
tually be the same singularity 1/r. It arises from the 
3D free space Green's function and becomes dominant 
as r — > away from the confining boundaries (see sec- 
tions IIV Bl and IVIII . This latter identification is one of 
the elements that will allow us to incorporate all par- 
tial waves I > 1 besides the s-wave. In addition, the 
above conclusions concerning CIR under parabolic con- 
finement will be seen to hold not only for general in- 
teractions V{r), as first examined in |3g|, but also for 
general non-parabolic confinement U{p) . In this sense 
the present formalism is of a more general character. 



IV. SCATTERING PHASE-SHIFTS 

In the following, a general formalism is developed in or- 
der to calculate the scattering amplitudes for any given 
finite short-ranged and spherically symmetric scattering 
potential V(r). This is done by expressing the amplitudes 
in terms of the scattering phase-shifts Si associated 
to this V{r). Our formalism is also able to deal with a 
broad range of confinements U{p) as well as scattering 
energies above the transversal ground state. Particularly 
relevant is the straightforward accounting of angular mo- 
menta couplings due to the confinement. We begin by 
providing the connection to scattering in free space. 



A. Phase shifts in free space scattering 

The standard formalism for free space scattering by 
a central potential is recast in a form suitable for later 



comparison. The main idea is to express all quantities 
relevant for scattering in terms of the phase-shifts. The 
rationale is that these phase-shifts can be considered as 
intrinsic to a given scattering potential V{r), without 
much regard to details of external boundary conditions. 

In the absence of any confinement, the scattering so- 
lution can be written as 



$(r) = $,(r) - J (fr' G{r, r') v{r') $(r'). 



(7a) 



where $i(r) is the incoming state and G{r, r') is the free 
space Green's function. In terms of outward and inward 
scattered waves, the latter reads 



G(r,r')=7+ 



Ak\r—r' I 



47r r — r 



^ — ik\r—r' \ 



(7b) 



for some constants 7± obeying 7+ + 7_ = 1. A single 
constant 7 can be introduced by setting 



7± 



1±7 



(7c) 



In free space scattering, the inward component propor- 
tional to 7_ is usually absent (i.e., 7=1). However, one 
expects this inflow of particles if there is e.g. an exterior 
confinement that forces the scattered particles back to- 
wards the center. This interpretation will be validated in 
section HV Bl (see also Fig|21l. $i(r") and G(r,r') are de- 
termined by the boundary conditions imposed on <&(t'), 
e.g. that of an axially symmetric solution. 

The scattering variables of interest here are the am- 
plitudes /out, in of outward and inward waves, scattered 
along a direction given by and cj) at large distances 
r ^ Rv from the scattering center. Using the expansion 
|r — r'l r — r ■ r' jr in Ea. (|7b|l . Ea. (|7a|l becomes 

$(r)«$,(r)-f /out— + /in^, r»i?y, (8a) 
r r 



with 



/o 



_7± 
47r 



i;(r')*(r'), (8b) 



where k = fcr/r depends on the direction (6*, 0) of r. It 
can be clearly seen that these amplitudes have a shape 
similar to their equivalent effective ID amplitudes /j^ in 
Eq.(jlE|l- To express the amplitudes in terms of the phase- 
shifts an angular momentum decomposition of /out. in is 
necessary. Following standard procedures (IH , Eg . we 
obtain from Ea. (|8b(l 

00 

/out.in = -7± 5I(T*)'(2^ + 1) Ti PiicosO) , (9a) 



1=0 



with the definition 



ri = ^J d^r' [jiikr') Pi{cos0')] v{r') $(r') , (9b) 
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where ji and Pi are the spherical Bessel functions and 
Legendre polynomials, respectively. The next step is to 
relate the elements t; to the phase-shifts. Decompose 
<i>i(r) and G(r, r') in spherical coordinates. Substituting 
then into Ea. (|7a|l for r > r' ^ Ry and carrying out the 
(/('-integration one obtains for r > Rv 



$(r) = Y.{2l + l){(3i-ik-fTi)ji{kr)Pi{cos9) 

1=0 

oo 

+ ^(2Z + !)(/?; + fcrO nz(fcr) (cos 0) , 
/=o 

having used Eq. (f7c|) and Eq. Ij9b|) and the constants (3 and 
P' are related to $i(r). On the other hand, the phase- 
shifts arise when directly decomposing this same solution 
<i>(r) in spherical coordinates and looking at the radial 
solutions outside the interaction region r > Ry, 



^{r) = [cos Si ji{kr) - sin (5; ni{kr)] Pi{cos6). 

1=0 

Here the phase-shifts Si are calculated from the behaviour 
of the radial parts Ri{r) inside the scattering region 
r < Ry only, irrespective of the details of the outer 
boundary conditions The boundary conditions are 
matched by properly choosing the constants q, i.e. by 
superposing the spherical components. Comparing these 
last two forms of <i>(r) yields. 



Tl 



(21 + 1)- 



cos Si — i"f sin Si 
1 Pi sin Si + P'l cos Si 
k cos Si — ij sin Si ' 



Z = 0,1, 



(10a) 
(10b) 



Considering the phase-shifts Si (or the scattering poten- 
tial v{r)) as given, the second relation Eg. (| 10 bp solves 
for the amplitudes /out, in hi Ea. (|9a|l (and also the inte- 
gral equation Ea. (|7a|l at distances r > Ry). The first 
relation Ea. (|10a|) in turn determines how to superpose 
the spherical components in $ = ciRiPi in order to 
match the boundary condition contained in Ea. H8a|) . 



B. Green's functions under confinement 

Just as the free space 3D scattering amplitudes in 
Eqs.® can be written in terms of the phase-shifts via 
Eqs.Q and Ea. (|10b|) . a very similar procedure can be 
applied to calculate the amplitudes characterizing the 
ID effective dynamics. In this section a key element of 
the formalism is worked out, namely, the Green's func- 
tion Gc under confinement. It is shown that Gc can be 
expressed in terms of the free space Green's function G 
as in Ea. (|12() . The final purpose, to be achieved in sec- 
tion IIV CI is to express it as a series in spherical coor- 
dinates, valid at least in the center of the guide where 
the spherical symmetry prevails over the cylindrical one 



due to the main assumption of a short range scattering 
potential v(r) and a confining potential u{p) that is suffi- 
ciently flat at the center, Ry ^ Such a series allows 
then the straightforward introduction of the phase-shifts. 

Consider the solution to Eq. (|4a|l given by the expan- 
sion in Eq.|0cj- 0^.0 can rewrite Eq.l^cJ in the integral 
form (keeping again the ^'-integration as in Ea. (j9b|) for 
convenience) 

*(r) = ^, {r) - J dPr'Gcir, r')w(r')*(r') , (11a) 



where the incident state is 



*^(r) 



TLE 



(lib) 



and Gc, the axially symmetric Green's function of the 
confining potential, has the form 

oc 

Gc(r, r') - ^rXpXip') Gn{z - z') . (11c) 

ri=0 

It is convenient to introduce a non-axially symmet- 
ric Green's function Gu under the confining potential 
u{p). This Gu then satisfies [V^ - u{p) -f fc^] G„(r, r') = 
—Sir — r'). If the (/(-independent solution \E'(r) is writ- 
ten as an integral equation using this (/)-dependent Gu, 
it follows from comparison with the integral equation 
Ea. (|lla|) that Gc and G„ are related by 2tt Gc{r,r') = 
J d(j)' Gu{i',r'). Compare now the differential equa- 
tions of Gu and that of the free space Green's func- 
tion G of Ea. (|7b|l in the region of the scattering po- 
tential such that r, r' <C Ru- Since u{p) « by as- 
sumption, they should differ at most by a homogeneous 
term, say A„, so that G„(r,r') sa G{r,r') + Au{r,r'), for 
r, r' <^ Rjj , where A„ satisfies the homogeneous equation 
(V^ + fc^)A„(r,r') = and we obtain 

Gc{r,r')^ J ^G{r,r')+Ac{r,r'), r,r'<^Ru, (12) 

where Ac = / dcj)' lS.u / 2tt . In order to obtain the function 
Ac and 7 in Ea. lffcjl in terms of the scattering param- 
eters, we first write the closed channel [n > ue) wave- 
functions of Gc in Eq. Ullcfl and G in Eq. Ij7b|l using cylin- 
drical Bessel functions. The latter represents a good ap- 
proximation to these transversal eigenstates for a given 
potential U{p) close to the guide center p <^ Ru where 
U{p) ~ 0. By subsequently comparing Gc and G we can 
determine Ac. 

First we express the transversal eigenstates (pn{p) h^ 
terms of the normalized Bessel function according to 



Nn 



TT^^RU 



Jo{qnP) , P < Ru, n > HE, (13a) 



where Nn = | Ji(r„+i)| ^ and r„+i is the (n -I- l)-th root 
of Jo, related to the transverse momenta (/„ by qnRu = 



7 



r„+i . A good expression for the roots r„+i is given by the 
cosine approximation Jo{x) ~ (2/7r.T)^/^ cos (x — 7r/4). 
This leads to 



3\ TT 



(13b) 



The corresponding energy dispersion relation e„ ~ is 
quadratic with respect to n compared to the linear one 
e„ ^ 2n for, e.g., 2D harmonic oscillators (see e.g. [ISl)- 
This should not cause any confusion since this excited 
spectrum n > is summed over in Eq. (|llc|) and leads to 
the confinement independent free space Green's function 
G in Ea. (|12|l . It is the lower part of the spectrum n < he, 
in turn, that relates to the term A^. 

Indeed, we substitute Ea. (|13a|) in the series for the 
closed channels in Eq. Hllc() yielding for r, r' <C Rjj 



G,(r,r')«£^„(p)^:(p') 



i-\yk^—q^\z — z' I 



1^ -^MqnP)MqnP')—^==^, 



7l — l-\-7lE 



where Sqn — t^/Ru is the increment of qn and we used 
|Ji(r„+i)p — 2/'KqnRu- For r,r' <C Ru a continuum 
approximation to the closed channel series can now be 
applied. Separating the real and imaginary parts of Gc 
(supposing fn{p) real for n < ue without much loss of 
generality) gives 



Gc{r,r') 



COS (^^/k^ ~ ql\z - z'\^ 



Tl = 



-^^n(.p) fn{p') 



sin (^^F - ql\z - z'\^ 



n=0 

qdq 



qi+n 



g-V q'^-k'^\z-z'\ 

, Mqp)Mqp') ==^,(14) 

47r „ 



where qi+UE > ^ designates the transversal momentum 
of the first "virtual" or closed channel n = 1 + ue- As 
for G, one uses now the identity , ^45.] 



^ik\r—r' | 

Airlr — r' 



^ (2 - (5o,m) cos [m{(p - (j)')] 



m=0 



qdq giy/k^^q^\z-z'\ 
X / - — J„i{qp)Jmiqp')- 







and the correct branch < Arg^ fc^ — (p- < tt in order to 



obtain for the integrated free space Green's function 
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k'^ — (p\z — z' 



qdq 

^Mqp)Mqp')— ^y^—, 



''qdq '^i'^ 
-^Mqp)Mqp') 



fc2 



q'^\z — z'\ 



-^Mqp)Mqp )—j==^. (15) 



To compare with the Green's function Gc in Ea. p4|l . we 
specialize Ea. (|15|) to r, r' <C Ru and impose the condition 
for low total scattering energy. 

The value of 7 follows from the comparison of the imag- 
inary parts of Eas. (|14() and l|15() . Changing the integra- 
tion variable from q to ^Jk"^ — q^ in the first integration 
on the rhs of Eq. H15() and using the low energy condition 
in the imaginary parts of Ea. (|14|l and Ea. (|15|l . 7 can be 
identified in the limit r, r' ^ 



7 = 27r 



t^k^Jk!^ -ql 



In other words, since r, r' <C Ru ^ l/k , both imaginary 
parts are approximated by a single and the same constant 
jk/An, independent of r and r'. Although exact values 
of |(^„(0)p can be employed for a given potential U{p), 
it is more instructive to use general approximations such 
as |</J„(0)p ~ l/nRfj obtained through normalization. 
An estimate that takes into account the spatial variation 
and the few nodes of is the square well approximation, 
namely, |<pn(0)p « N^/ttR1j. Therefore, without loss of 
generality, one can thus write 



E 



2N' 



^^Rlk^/W^l 



(16a) 



It should be mentioned that one can improve on how 
these low energy poles g„ of Gc are treated. However, it 
suffices to use for now the above result in order to present 
the formalism. 

By comparing the real parts of Eq. 1)14(1 and Eq. H15() , it 
is seen that Ac in Eg. 1(12(1 is given by 

Ae(r,r') 



91+11 



qdq e-V^I^-'l 

-T—Mqp)Mqp ) , „ ,„ 

47r _ /j2 

^Jo{qp)MQP') sm{p\z - z'l) 
47r 



^fnip) ^n{p' 



, sin (fc„|z - z'l 



2k n 



8 



where fc„ = yjk'^ — and q = yjkp' — . The first term 
on the rhs stems from the offset between the lower lim- 
its of integration qi+UE ^^^^ ^ iii Ea. (|14|l and Ea. (|15|l . 
respectively. It accomits for the discreteness of the low 
lying transversal states due to the confinement. The sec- 
ond term in the square brackets is of order (fc^/87r — 
^^^Q A^^/27ri?^) |z — z'\ and thus is relatively smaller 
by a factor of |z — z'\lRu compared to the first term, 
whose order of magnitude {q^j^^^ — 
can be estimated by calculating the integral for r, r' . 
Neglecting this second term yields 



A,(r,r') 



^ jQ{qp)Jo{qp')e' 



(16b) 



where q = (fc^ +p2)i/2 ^^d Pc = (g?+„E ~ k'^)^/'^. The 
relation Ea. (|12l) is then proved and will be the basis of 
an angular momenta decomposition of the effective am- 
plitudes in section HV CI In this regard, the weak de- 
pendence of Ac on the coordinates r and r' is kept, since 
this dependence introduces couplings between orbital an- 
gular momenta. This will be seen as Ac is decomposed 
in spherical coordinates in the following. 



C. Angular momenta decomposition 

We are now in the position to express the effective 
ID scattering amplitudes in terms of the scattering 
phase-shifts Si of a given scattering potential. The first 
step is to decompose these amplitudes in Ea. HSbl) in a 
similar way as was done in the free space case in Eqs.Q 
by introducing "matrix" elements similar to r/. Subse- 
quently these quantities are related to the phase-shifts di 
in a manner analogous to Eas. (|10|) . 

Applying the condition Ry <C Rjj, the transversal 
states (fin in the integrand of the amplitudes Ea. lj5b|) 
should be well approximated by the Bessel functions 
Ea. (|13a|) . This is because r' ^ Ry whereas |fc„| ~ 
Qn ~ ^/Ru- It is enough then to find expansions of 
e^"^^"^ Jo{qnp) in spherical coordinates, with = q^+^ 
For this purpose one may invert the following identity |4 



then express the amplitudes in Eq. Ij5b|l as the follow- 
ing series in angular momenta for n < ue 

ft - fug ± fnu (18a) 

47r {21 + l)ain ^ 47r {21 + 1) ain 



E 



ti±y: 



I even I odd 



T, 



having used the parity property Pi{—x) = {—yPi{x). 
The "matrix" element T; is defined by 



T, 



1 



V An 



d-^r' [ji{kr') PiicosO')] v{r') l'(r'). (18b) 



The non-zero momenta / > 1 contributions to ft stem 
from the (small) dependence on z' and p' of the integrand 
in the definition of ft ■ Their neglect amounts to assum- 
ing a point like zero range interaction v{r), for which only 
the s-wave remains. In the free space scattering case, the 
corresponding quantity equivalent to this T/ is r/ given 
in Ea. (|9b|l . The key difference now is that each T/ will 
depend on all other Tis, meaning that the confined scat- 
tering solution ^(r) in the definition of T/ couples the 
angular momenta, while the free space solution <I>(r) in 
the definition of r; does not. 

The next step is to relate the T; to the phase-shifts 
5i. From the results in Eas. (|16|l for the approximation 
Ea. (|12|) . the Green's function Gc can be decomposed into 
spherical coordinates in the region r',r <^ Rjj by sepa- 
rately decomposing G and Ac. The free space part G can 
be decomposed following the procedures of section ITVAI 
Then, for r' < r. 



— G(r,r') = ik'^ji{kr')['yji{kr) 

1=0 



i ni{kr)] 



21 + 1 

x—^Pi{cose)Piicoae'), (19) 
47r 

having used the expansion of £±''=1''-'' l/47r|r — r'\ in 
spherical harmonics and 7± — (l±7)/2. To evaluate 
Ac we use the analytic continuation 9 ^ tt/2 — i9' of the 
first identity of this section to obtain 



Pi{cos0)ji{kr) 



2*' Jo 



ikr cos 9 cos w 



dw sin w e 



X Jo{kr sin 6 sin w)Pi (cos w) 



by using the orthonormality property of the Legendre 
polynomials. As a result, one obtains for r <C Ru 

oc 

e*^-"^y.„(p) = ^z'(2/ + l)az„ji(fcr)P,(cos^), (17) 

1=0 

where ~ NnPi{kn/k)/TT^^'^Rij. For more quantita- 
tive results, these coefficients a;„ can be evaluated nu- 
merically from the exact eigenfunctions ipn{p). One can 



e-PVo(gp) = + l)mp/k)\ j i{kr) Pi {cos 6) , 



1=0 



where p ~ A;sinh6'' and q — k cosh 6' such that k^ — 
q^ —p^. The expansion of Ac in Ea. ljl6b(l for r, r' ^ Rjj 
can then be written as 



Ac(r,r') = -f; p ^Pi{a,,np/k)e'^^^'^^' Mqp') 



1=0 



X i'{2l + 1) ji{kr) Pi {cos 9), (20) 



where azz' = sign(z — z'). Inserting Ea. (|2()|l and 
Eg. lfT^ back into Ea. ltT^ . the scattering solution ^'(r) 
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in Ea. ljllall becomes for Rv ^ r ^ Rjj 

oo 

*(r) « (2/ + 1) [a, + 7, (z)-i7A:T,] J, (fcr)P, (cos (?) 

oo 

+ ^ ^' (2? + 1) [ A:Ti ] m {kr)Pi (cos (?) , (21) 
;=o 

where ai = X]"=o ^nOLin comes from the expansion of 
e^^^^ ipn{p) in the initial state ^'^(t) in Ea. ljllb|) and7;(z) 
is defined by 

e-""'-' Jo(gp')] v{r')^>{r'). 

The z-dependence of this ^i{z) arises from the sign of 
Gzz' in Ac. Because of this z-dependence, the above de- 
composition of ^'(r) is not fully spherically symmetric. 
It reflects the fact that an approximate radial solution on 
the left side of the guide at some zi — — |zo| needs not co- 
incide necessarily with one on the right side at Z2 — \zo\ 
since the problem has no perfect spherical symmetry. 

A better understanding of the role of this 7i(z) arises 
when one tries to write it in terms of the T;'s that de- 
termine the amplitudes in Ea. (|18a|) and the solution 
in Ea. (|21|l . Expanding e°'"'P^ Jo{qp') (see above) in the 
integrand in the definition of 7i(z) one obtains 



7i(^) = + 

47r 



dpPi{ip/k)P,{ip/k) 



'js{kr')P,{cos9')]v{r')^{r'), 



assuming the summation and integration can be in- 
terchanged (see section IV C(l . If it happened that 
(czz')'^^ = li then 7i(z) would be a linear sum of the 
elements Tg and the scattering solution in Eq. (|21() would 
have spherical symmetry. Such a constant 7i(z) could be 
possible if one could neglect in the series for 7i(z) terms 
for which the parity of Pg is different from that of Pi, 
such that for each I, only terms satisfying I + s = even 
would be retained. This condition can be related to 
the fact that the confining potential u{p) only couples 
angular momenta with the same parity, in the sense 
that {l\u\s) = J de'Piicos9')u{rsm9')Pg{cose') if 
I + s = odd. Therefore, a consistent approximation to 
the scattering solution \['(r) is found if one neglects odd 
parity combinations in the series defining 7z(z) and sets 
7;(z) « 7;, such that 
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s[l] 

= E(2s + 1)P,,T, 

s[l] 



Vc/k 

k I dx Pi{ix)Ps{ix) 





T, 



0,1,2,..., (22) 



where s[l] denotes the sum over even (odd) s for even 
(odd) I. 

Using Ea. (|22|l . the scattering phase-shifts 5i can now 
be introduced, since the expansion in Ea. (|21|l becomes 
spherically symmetric. Indeed, this expansion can be 
conveniently rewritten by the introduction of constants 
c'l and 5i defined, for the moment, by (/ = 0, 1, 2, . . . ) 



c'lCOsSi = i\2l + l)[ai+"ii-i"fkTi\ 
c'lSinSi = -i\2l + l)kTi, 



(23a) 
(23b) 



such that our solution has the form 



oc 

c'l [cos 5i ji (kr) - sin 5i ni (fcr)] Pi (cos 0) . (24) 

1=0 



In order to identify the Si, one notes, on one hand, that 
this wavefunction ^'(r) can be directly obtained as a 
spherically symmetric series c;i?;(r)P; (cos 6*) in the 
region r Rjj, where u{p) is negligible, if one solves 
the Schrodinger equation moving outward from the ori- 
gin r = (see Fig 15}. The radial part Ri can then be 
separately determined and for Ry r Ru it should 
be given by Ea. l24|l . provided one chooses q = c\. On 
the other hand, the solution for the scattering potential 
v{r) without confinement is given by a series with the 
same radial parts Ri but different constants ci ^ c'l (see 
Section FlV A|l . The asymptotics of Ri in this free space 
case is then also given by Eq. H24|l so that 5i are indeed the 
phase-shifts of the unconfined scattering problem. The 
difference between the confined and free space solutions, 
in the region Rv ^ r <C Ru, is thus accounted for by 
the different constants cj and ci, respectively, related to 
distinct boundary conditions. Besides, 5i depends solely 
on the solution of Ri in the interior region r < Ry 463. 
The above defining relations for C; and 5i in Eas. l|23(l are 
then equations that determine c\ and T; in terms of 5i, 
namely. 



Cl 



(2/-f l)(a^+7^)^' 



Ti 



ai +7/ 



cos 5i — sin Si i^k — k cot Si 

Finally, using Ea. (|22|l one obtains the matrix equation 
relating Ti to the full ensemble of 3D free space scattering 
phase-shifts Si (Z = 0, 1, 2, . . . ) 



(i-fk - k cot Si) Ti^ai + V(2s + l)PisT, 



s[l] 



(25) 



The coupling of angular momenta brought about by the 
confinement is a result of Ac, which accounts for the 
confining geometry and the discreteness of the low lying 
transversal states that should be resolved at low energies. 
The main equations that allow us the analysis of the ef- 
fective ID scattering in confined geometries are Eqs.(|SJl, 
Ea. ((TKal) and Eq.^. The probabihty conservation con- 
dition Eq. (|5c|) serves to gauge the range of validity of the 
results. Although the above result is enough to present 
the formalism, improvements to Ea. H25|) can be system- 
atically made if necessary, e.g. by better dealing with the 
poles Qn in Ea. ljltiall and the constants a;. 
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V. CONFINEMENT INDUCED EFFECTS 

We now use the above formalism to analyze confine- 
ment induced phenomena that occur both in the scat- 
tering by a central field and in collision processes. Al- 
though related to each other, two important phenomena 
can be distinguished: resonant scattering (sections IV HI 
and lVDl and weakly localized states induced by the con- 
finement fsection fVEl) . In the present paper, our formal- 
ism is focused on the single mode regime. It is shown 
also that previous results of atom-atom collisions under 
parabolic confinement are recovered, re gard ing not only 
the scattering resonances (CIR) [SJ, [Sa, |3g but also the 
confinement induced weakly bound states (s^. CIR in 
atomic collisions involving excited transversal states at 
energies > qf are treated in l35l for interactions in 
the s-wave approximation and in pq for general V{r) . 
See also independent related works for 2DEG mesoscopic 
systems [H H El El . 



A. Single mode regime 

In this regime, the total scattering energy E = h^k^ 
allows for only the transversal ground state (poip) to be 
effectively occupied k'^ = qQ + k^ < qf, i.e., to be an open 
channel. All excited states can only be virtually popu- 
lated due to the scattering of the incoming state. Then 
He = and 6„ = dn,o- In terms of the even fog and 
odd /oti angular momenta decomposition of given in 
Ea. (|18al) . namely Jq = fog ± /ou, the current conserva- 
tion can be written as 

(Re{/oJ + l/ogi') + (Rc{/on} + |/o«|') =0. (26) 

A suitable parametrization of /og.u is obtained by intro- 
ducing ID scattering parameters 6g and Su, such that 



fog 

fou 



1 



1 + i cot Sg 

1 

1 + i cot Su 



(27a) 
(27b) 



In this form, Ea. H2t)|) is verified if 6g^u are arbitrary real 
numbers, since then Re{/og,«} + |/o£(,Mp = vanishes sep- 
arately. There is by now no a priori reason for the param- 
eters Sg^u to be both real. However, in the case of general 
potentials v{r) and u{p) considered in this work, it will 
be seen in the next sections that Eas. (|27|) with real Sg^u 
provide a valuable means to establish the boson-fermion 
and fermion-boson mappings as well as the conservation 
condition Ea. (PS|l . 



B. Boson-fermion and fermion-boson mappings 

Consider now a coUisional process between identical 
particles. It is convenient to split ^>{r) into symmetric 



'ipg{z)ipo{p) and antisymmetric ipu{z)ipo{p) parts, with re- 
spect to r ^ — r, 



*(r) = 



Jkoz 



± ikQ\z\ 



foip) 



[ipg{z) + ipuiz)] tpaip) 



\z\ > R[, 



ipg{z) and ipu{z) can be deduced by using the form 



Ea. (|18a|) for the amplitude /q in terms of its even fog 
and odd fou sectors. For the symmetric part, we obtain. 



V'g(^) = (l + /og)cos(A;o2:) + «/ogSin(fco|z|) 

= e'^^cosikolzl + Sg), \z\>R[^, (28a) 

where Ea. (|27al) has been employed. The antisymmetric 
part is given by 

tpu{z) = i(l + /o„)sin(fco2;)±/o„cos(A:o2;) 

= ie'^^ sin {koz ± 5^) , \z\ > , (28b) 

where the plus (minus) sign in ^puiz) refers to z > 
{z < 0) and Ea. (j27b|l has been used. 

If the colliding particles are e.g. spin polarized ^'(r) 
must be (anti-) symmetrized. For bosons, the correct 
effective ID scattered wavefunction in the asymptotic 
region \z\ > R[jj is ipgiz). The resonance condition 
can be identified as fog = —1 or Sg ~ and it 

follows from Ea. (|28a|l that the wavefunction becomes 
t\}g(z) = — zsin (fco|z|), which is with respect to its modu- 
lus the wavefunction of a pair of non-interacting fermions, 
as can be seen by setting /o„ = (or (5„ = 0) in Ea. (|28b|l . 
This is the well-known fermionic mappi ng o f strongly in- 
teracting, inpenetrable bosons in ID The 
bosons would not be allowed to be located at the same po- 
sition 2 = 0, supposing the asymptotic solution Ea. (|28a|) 
could be extended towards the origin Izj <C R!^^ ■ 

For fermions, a reciprocal mapping exists at resonance 



fo-u 



k/2| 



(29) 



From Ea. lj28b|) . the wavefunction becomes tpuiz) — 
=Fcos(fco2;), which is the wavefunction of a pair of non- 
interacting bosons [3^ . Clearly, from the exact, fully 
antisymmetric wavefunction of a pair of spin polarized 
fermions, the probability of finding both of them at the 
origin r = is strictly zero. The bosonization implied 
by this confinement induced resonance means that the 
asymptotic behavior of the fermions far from the inter- 
action region is that of free bosons. 

These mappings between identical colliding particles 
are valid also under a longitudinal confinement along the 
z-axis. In order to neglect couplings to the center of 
mass, this longitudinal confinement can be chosen to be 
parabolic. The characteristic oscillator length should 
then be flat enough, such that ^ R'id- The asymp- 
totics of the symmetric and antisymmetric eigenstates 
with (positive) energies Eg,u = ^^^g m/^M are given by 
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Ea. (|28a() and Ea. lj28b|) . respectively. The correspond- 
ing eigenvalues fc^ „ = 9o + ^Og,u ^e estimated by 
solving the equations '03,«(-R||) ~ for fcog.u, namely, 
cos(fcogi?|| + Sg) = and sin(fcoMi?[| + (5„) = 0. If the 
parameter Sg = 5g{kog) is made equal to ±7r/2 by vary- 
ing, e.g., Rvi Ru, or one sees that ipg{z) is once 
more mapped onto an antisymmetric eigenstate with- 
out interaction v{r) = (or Su = 0), as already seen 
above for free longitudinal motion. Analogously, i^u^z) 
is mapped onto a non-interacting symmetric eigenstate 
as 6u = ^«(fco«) ~^ ±7r/2 . This latter type of bosonic 
mapping of fermionic eigenstates under longitudinal con- 
finement was first verified numerically in '301 for the col- 
lision interaction V{r) ~ d/ cosh(r/6) . 

Note that, if one is not restricted to the case of colli- 
sions, this result can be formally quite general regarding 
not only the interaction potential v{r) , but also some suf- 
ficiently flat transversal and longitudinal confinements. 



C. The effective amplitude /q 

The considerations in the preceding section assumed 
the general form Eas. (|27() for /og,M . In this section, 
is explicitly calculated by solving the matrix equation 
Eq.(|nnil for the elements Ti . From k ^ 1/Ru, it fol- 
lows that kRv ~ Rv/Ru ^ 1, which is the condition of 
low scattering energy in 3D free space scattering and the 
phase-shifts 



tan (5; = tan(5;(fc) ~ k 



21+1 



1/Rlj+' 



are generally small 46] for large Rjj (or small Rv)- One 
expects then that orbital angular momenta higher than 
the leading contributions, e.g. the s- and p- waves, should 
not significantly change the main features arising from 
the latter. Indeed, the calculation of fog and /ou can 
be done separately, since the even and odd angular mo- 
menta in Ea. (|25l) are uncoupled from each other. In the 
single mode regime, Ea. HlSal) becomes 7 = 2/df, k ko and 
du = Ru/No, where fco = {k'^-QoV^'^ and in Ea. ((TSat 
depends only on the ratio ti = Ti/k^. The matrix equa- 
tion Ea. (|25|l in this case can be rewritten in terms of 
these tl . Thus for all / , 



{2i- [kcoiSi + {21 

= aidfj 



l)Pii]dfjko}ti 



kodf 



(2s + l)Pu t. 



where the diagonal term s[l] — I is excluded from the 
series on the rhs. Assume now in this equation for 
tl that this series converges. As I increases, ai — 
aio = Piiko/k)/n^/^du (see Eq.^ and Eq.^) de- 
crases, while Pis (see Eq.(|521)) should decrease mildly 
due to the oscillation of P; (ix) . Then the rhs of this 
equation for ti should decrease relatively smoothly with 
I . On the other hand, the phase Si decreases exponen- 
tially Si ~ l/R^}^^ with I at low energies, so that cot Si 



increases exponentially. Then ti on the Ihs should be ex- 
ponentially small with increasing / . This is consistent 
with the assumption of convergence of the series on the 
rhs [see also discussion on 7; (2) after Ea. H21|l ]. Apart 
from this, one can use the smallness of the phase-shifts 
1^2! ^ |<5o| and l^s] <C \Si\ in order to explicitly diagonal- 
ize the finite subspaces I = 0,2 and I = 1,3 separately. 
Then the couplings to I = 2 and I = 3 are seen to have 
a minor effect to the leading terms I — and / = 1, 
respectively. Hence we put 



tl 



jdjjai 



[kcotSi + {2l + l)Pii]ldlko 



l = Q,l,. 



The even sector /og is then given by the leading term, 
e.g, the s-wave, for which fccot(5o ~ — 1/a, where a is 
the s-wave 3D free space scattering length From 
Ea. (|18a|) one obtains 



Og 



1=0.2,. 



4tt {21 + 1) a;o 
2i 

1 



1 + i 



2a 



aPni 



ko 



(30a) 



The odd sector /qm is approximated by e.g. the p-wave. 
In this case, the relevant quantity j47| is the scattering 
volume Vp related to Si by k cot Si = — 1/V^fc^ . Then 



fov 



^ An {21 + 1) aio ^ 
1^ 



i=1.3, 



2i 



3[Pi(fco/fc)]' 



1 



kn 



(30b) 



The range of validity of these results can be obtained 
by imposing e.g. the conservation of probability. As 
discussed in connection to Ea. H26() and Eas. H27() . the ID 
scattering phase-shift Sg is real for the result Ea. H3Ua|l . 
Then a sufhcient condition for Ea. (|26|l to hold is 

\h)\<^k^l/Ru (31a) 

since then 3[Pi(fco/A:)]^ ^ 1 so that /ou ~ and the 
parameter Su ~ can be considered also real. In this 
first region of low longitudinal momenta, the angular mo- 
menta couplings in the matrix equation for ti should be- 
come negligible due to the factor fco that multiplies the 
series X]s[/]^;(^'^ ^)Pists containing other angular 
momenta. The present uncoupled solution to ti should 
then be a good approximation. The second region is at 
relatively large momenta, when 2,[Pi{ko/k)Y ~ 1, i.e., 
when 



M^qol^^l/Rl 



(31b) 



since then S^ can also be taken as real. The decoupling 
between odd angular momenta rests then on the small- 
ness of the phase-shifts Si . These are, nevertheless, the 
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most interesting regions for the momentum fcg that are 
relevant to resonant scattering in symmetric and anti- 
symmetric states, respectively. For more detailed appli- 
cations these regions can be enlarged by improving on 
Eq.lESl). 



D. Confinement induced resonances - CIR 

As presented in sections lV Al and IV Ijl in general terms, 
the CIR is one remarkable effect arising due to the pres- 
ence of a strong confinement. In the following, it is ex- 
plicitly expressed in terms of the parameters determining 
a given scattering potential and a given confinement. 



1. Even angular momenta 

The symmetric wavefunction in the single mode regime 
is 'ipg{z)(pQ{p) and only the symmetric sector /og con- 
tributes. It can be used to describe e.g. spin-polarized 
bosons or fermions in an antisymmetric spin state. In the 
expression Eq.JSOijl for this /og, Pqo is given in Ea. l|22|l . 
thus 



rPa/k 

k / dx PQ{ix)PQ{ix) 
Jo 



(32) 



having used the definition of pc in Eq. (|16b|l and the single 
mode condition — . 

At small scattering lengths |a| ^ Ru , one can neglect 
this factor Pqo in Ea. (|30a|) . since then |a|Poo ^ 1 , and 



1 



Og 



1 



4 

2a 



\a\ < Ru 



This is also true at "high" energies | fco I when k qi such 
that Poo ~ 0. One must, however, be careful about the 
validity of the present approximations in this region of 
large momenta |fco|, as mentioned above in the context 
of Ea. lj31b|l . The only possibility for resonance /og = — 1 
is then at zero energy fco = . 

At large scattering lengths Ry '^\a\ ^ Ru , this factor 
Poo plays an important role in the scattering process at 
low momenta |fco| ^ ^/Ru for which Pqo ~ y/ — qo is 
maximum. Indeed, one obtains then, for |fco| <C i/Ru i 



fog 



1 



2a 



1 



C'g 
du 



(33a) 



where the constant C depends only on the first closed 
and the last open channels, namely 



c 



du\/qf-qo 



Y7^ Ve]~~eo- 



(33b) 



Recalling that /q k, fQg, since fou ~ for such small 
momenta |fco| (see Eas. (|3U|l '). this result implies that the 
scattering can be described by means of an effective ID 
interaction potential Vid{z) of zero range 



Vid{z) = giD^iz) . 



2a 



1 



giD 



C'a/di 



In other words, Ea. H33a|l is the solution to for the hy- 
pothetical pure ID scattering problem under this poten- 
tial ^113(2:) and is valid also for distinguishable particles. 

As discussed in section IIIII for atom- atom collisions, 
this same ID potential Vioiz) at low momenta |fco| was 
first shown in to be a direct result of parabolic con- 
finement and of the 3D zero range interaction potential 
Vs{r) of Ea. l|t)a|l . The occurence of this coUisional reso- 
nance at a± = Ca for general interactions V{r) was then 
shown in 3^] . 

For other types of scattering processes, e.g. in meso- 
scopic 2DEG systems [H El [S [H, Ea. (|SSal) shows 
therefore that at low longitudinal momenta fco the exis- 
tence of a resonance of the type Jq ~ fog = — 1 is char- 
acterized by an infinite effective ID coupling strength 



\9id\ 



du ~ C'a , 



(34) 



and is indeed not restricted to atom-atom collisions un- 
der parabolic confinement or to two-dimensional quan- 
tum point-contacts with a single impurity p^sf . It is a 
general effect of single mode scattering at low velocities 
under quasi- ID confinement. 



2. Odd angular momenta 

For odd angular momenta, the procedure is similar to 
the previous analysis. The antisymmetric part ipu{u) 
shows a scattering resonance when the associated am- 
plitude fou = — 1- At this resonance, the system can 
be mapped to a non-interacting symmetric wavefunction 
as discussed in connection to Ea. H29|) . The term Pn is 
given in Eq. H22|) , from which the actual condition of CIR 
in the amplitude /ou in Ea. lj30b|) can be worked out for 
given scattering and confining potentials. Note that due 
to the factor 3[Pi{ko / k)]'^ in the numerator of Ea. (|30bp . 
the odd component /ou is non-negligible only at relatively 
large longitudinal momenta |fco| ~ 1/Ru, in contrast to 
the even amplitude /og discussed above. 

For atom-atom collisions under parabolic confinement, 
the existence of CIR beyond the s-wave was first pre- 
dicted by neglecting couplings to orbital angular mo- 
menta and isolating a single partial wave I > and its 
phase-shift 61 as the sole dominant contribution [3^ . A 
direct consequence of contributions from non-zero partial 
waves Z > is, among others, the possibility to define ef- 
fective p-wave zero-range scattering potentials analogous 
to Vs{r) and Vioiz), but which act only on antisymmet- 
ric wave functions of a colliding pair [la . |49| . 
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E. Weakly bound states 

Besides CIR, a second effect of strong cylindrical con- 
finement is the prediction of weakly bound and localized 
states. For atom-atom collisions this implies the forma- 
tion of quasi- molecules. As first shown in jSTj for zero- 
range atom-atom interactions under parabolic confine- 
ment, these bound states exist even when no bound state 
occurs in free space. They have recently been observed in 
optical traps loaded with ^"K atoms '31^. Similar bound 
states localized around impurities are found in indepen- 
dent studies of mesoscopic systems 0, 0, 0, . In 
this section, a localized state for general cylindrical con- 
finement is calculated from the even, more precisely zero, 
angular momentum sector. 



1. s-wave binding energy and wavefunction 

For the free space I = angular momentum, the bound 
state of the (attractive) scattering potential V{r) has an 
energy Eb/ given by 



2 
1 

-1 

-2 
-3 
-4 
-5 
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Square-Well 

Free-Space 





-2 



la 



FIG. 3: (Color online) s-wave binding energy Eb for 
parabolic (dashed curve) and square well (upper solid curve) 
confinement, scaled by the respective ground state energies 
eq. The free space energy (lower solid curve) Esf is scaled 
with the square well value of eo and is defined only for posi- 
tive a. The dotted line marks the continuum threshold under 
confinement. 



Ea. (|^ for Poo, the equation for xb becomes 



EBf = -- 



2^ 



a:^ Rv 



whose relationship to the scattering length a via KBf ~ 
1/a holds only when a ^ Ry > [4^. At negative 
scattering lengths a < 0, one can speak of a "virtual" 
state close to be incorporated to the spectrum of V{r). 
Under lateral confinement and along the radial direction 
of the cylindrical trap, the tail of its wavefunction e^'^^J"' 
is changed to vanish at a finite distance from the origin, 
i.e. at the edge r — p — Ru of the guide when \a\ ~ Ru 
(see Fig. 0. This lateral squeeze lifts EBf < by an 
amount ep. It can be sufficient for this state to pass the 
threshold without confinement E = as Rjj decreases 
further. 

This confined localized state with energy Eb satisfies 
Eg . H4a|l with replaced by 2/ii?s/fi,^. This replacement 
is equivalent to redefining ko to be the imaginary num- 
ber koB = ± j('?o - 2pEB/fi^y/^. Since the wave e*'"'^'^ 
should be absent from the solution 'J and e**^"^!^' should 
show a bound state like exponential decay, l//(^(fcos) 
must vanish at a positive imaginary value of fcos , such 
that the interaction part e*'^°^l^l outweights e**^"^^. Let 
then koB = ixB/a, so that xb > if a > and < 
if a < . In terms of this xb , Eb can be rewritten as 



Ef 



2fi 



axB > ■ 



(35) 



From Eas. (|3U|l . a root oil/ arises either from the even 
sector 1/fog or from the odd sector 1//om. In order to 
recover the I = free space state Eb/, this root should 
come from the even sector Eq. (|30a|) . Using then the sub- 



stitution fc^ 2pEB/fT? 



^jar in the expression 



2s' 



Xb 



.s = — , (36a) 



where \a\ arises as 1/a^ is taken out of the square root in 
Poo ■ Another enlarged form is the quartic equation 



(C 



,2^2 



Wxi 



is^XB 



4s^ 



0. 



(36b) 



For each \s\, one should pick up the correct positive (neg- 
ative) root a;^ if s > (s < 0), as indicated in Ea. H35|l . 

The weak confinement limit is |s| — > 0. A finite root 
xb ~ +1 exists only for a > as can be seen from 
Ea. (|36al) . This means that Eb ~ Eb/, when Rjj oo 
or when a ^ 0+ as expected. For a < , a trial Taylor 
expansion xb ~ xis + X2s'^ in Ea. l)36bll yields Xi = and 
X2 = —2. A localized state with energy 



EB^eo-4: 



\E 



Bf\ 



\a\ du 



is thus expected also for negative scattering lengths. 
In this case, the confinement has been able to turn 
the "virtual" state EBf of V{r) into a real bound state 
localized around the scattering center. 

The strong confinement limit is \s\ ^ 1 (assuming the 
basic requirement Ry <C Rjj is not violated). Setting 
then Xb ~ Xoo/s~^ in Ea. (j36b|l and reexpressing it in 
terms of gives for positive and negative a , respec- 
tively, Xoo = ±{[(16 + C"4)i/2 _ C"2]/2}i/2. For both 
positive and negative roots xb , the same fraction of the 
transversal ground state energy eo is reached as |s| — > oo 



Ef 



\qodi 



eo : 



\a\ > dr 
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For parabolic guide du = a± [see Eq. (|16a|l with the exact 
value |(^„(0)p = I/tto^]. Then C' — 2 and qoa± — 
\/2, both obtained by directly using the exact energies 
£0 = ^^9o/2/^ = ftLO± and ei = = 3hLU± in 

terms of uj±. Hence Eb = (2 — V2)en = 0.586eo in good 
quantitative agreement with |3ll l37| . For square well 
confinement, it follows from Ea. ljl3b|) and Ea. lj33b|) that 
C" = (20/3)1/2 = 2.582 and qodu = (3/2)i/2 = 1.225, 
thus Eb — 0.631eo ■ These energies are plotted in Fig.|3| 
as function of the confinement length scale du- 

The spatial shape ^'_B(r) of this confinement induced 
weakly bound state in the region far from the scattering 
center is 



2. Higher angular momenta 

Localized states from higher angular momenta contri- 
butions and their couplings among each other can in prin- 
ciple be systematically obtained by solving Ea. (|25|l and 
analyzing the poles of the full amplitude — fog ± Jq^ ■ 
For instance, a d-wave pole should come from retaining 
the I = 2 contribution. This and other poles from the 
odd sector /o„ are treated in detail elsewhere. 



VI. DISCUSSION 



-XB\z\/a 



\z\ > R[ 



D ' 



(37a) 



setting aside overall constants and where is the proper 
root of Ea. (|36|l . This asymptotic form is common to both 
positive and negative scattering lengths a, with uxb > 0. 
Closer to the origin, the wavefunction can be calculated 
by analytic continuation of the intermediate asymptotics 
in Eq.lEU and Eq.(|22I)- There the function To = ro(fco) 
diverges at ko — ixb/o,, which is the condition for the 
pole of f^. After this analytic continuation, the step 
towards Ea. H24|) can no longer be taken (it was taken 
in order to calculate Tq — To{ka) for real ko). In this 
series Ea (|21|l , this divergence then singles out terms pro- 
portional to the constant To{ixB/a,) — > cxd . As a result, 
apart from overall constants, the wavefunction ^'^(r) has 
the following form (for Ry <C r ^ Ru) 



DC 



2/dl 



(37b) 

joi^Br) + Kb no{KBr), 



where = 2^EB/fi^ and Eb is given in Ea. (|^ . For 
Eb < 0, it can be seen that both imaginary values of 
Kb yield the same result. In particular, the free space 
bound state is recovered in the limit a/du 0+. Using 
Eb ^ Ebj ~ -T? I2[ic? in Ea. (|37b|l . one has 



^'B(r) cx 



a/di 



For negative scattering length a/djj 0~, the free 
space "virtual" bound state turns into a real one whose 
energy was calculated in the previous section, Eb ~ 
eo — A{a/du)'^\EBf\- The wavefunction Ea. (|37b|l becomes 
instead 



^B{r)oc jo {qor) + qo^o {qor) , a/ du 
a 



0" 



On the opposite limit of strong confinement \a\/du — > oo, 
Eb tends to a positive fraction A of eoi and the wavefunc- 
tion becomes a superposition of both terms jaiV^qaf) 
and no{^/XqQr). 



Here the role of the continuum approximation in sec- 
tion llV Bl is discussed by comparing with previous models 
for scattering under confinement. In the interior of the 
guide r, r' <ti Ru where the confining potential is negligi- 
ble, the Green's function Gc{r, r') under confinement ap- 
proaches the 3D free space Green's function G{r,r'), as 
shown in Ea. (|12|) . Therefore, the singularity of Gc{r,r') 
as r,r' ^ is not only coincident with (see section 
but it is essentially the well-known singularity l/\r — r'\ 
of the 3D scattering scenario without confinement. In ad- 
dition, one expects that the details of the confinement, 
whether parabolic or not, should not be important for the 
physical understanding. In fact, the spectrum of excited 
transversal states sums up to yield G{r,r') (integrated 
over the axial angle 0'), as the continuum limit shows. 
The only condition is a short-ranged scattering potential 
Rv ^ Ru or a correspondingly "flat" confinement. 

The remarkable effects induced by the confinement are 
in turn directly related to the low lying transversal states 
as also demonstrated in for a two-dimensional prob- 
lem. The discreteness of these states is then captured in 
the factor 7 and in the correction term Ac(r,r'), which 
introduce renormalizations and couplings of orbital an- 
gular momenta not accounted for by G(r, r') alone. 

For a more quantitative analysis, consider the specific 
case of atom-atom collisions in the low energy s-wave 
app roximation and under parabolic confinement [33 . l35l 
l37| . As discussed in section iHll the discrete summation 
over the transversal states can be dealt with in order to 
extract the singularity of the scattering solution and to 
determine the value of the constant C = 1.4603 ... in the 
effective scattering amplitude . On the other hand, 
the corresponding value of C in the continuum limit for 
this case is C = liuis^ooifQ ds' / ^/s' - X]s'=i W^') ~ 
Jods'/y^s' = 2 (see Eq.(9) of "34]), in agreement with 
G' = 2 calculated in section IVeI This value C" = 2 
arises by using the exact harmonic oscillator energies 
eo = hLUj_ and ei — 3hu!± . It can be improved on by 
e.g. numerically computing Ac. As already pointed out 

2n is 



in section IIV Bl the exact excited spectrum e„ 
substantially different from the "flat" potential approx- 
imation given in Ea. (jl3b|l . However, the discrete sum 
over this exact parabolic spectrum can be shown |3 [s^ 
to have an equal singular behavior l/\z\ as in free space 
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(see section Hence also within the discrete ap- 

proach, this demonstrates the major role played by the 
lowest transversal states, whereas the sum over the ex- 
cited states turns out to be qualitatively confinement- 
independent. 

VII. CONCLUSIONS AND PERSPECTIVES 

An analytical treatment of quasi-lD quantum scatter- 
ing by spherically symmetric and short but finite range 
potentials in general cylindrical confinement is devel- 
oped. The full scattering wavefunction is calculated non- 
perturbatively without partially resuming perturbative 
series. All phase-shifts of the scattering potential can be 
readily incorporated. This formalism provides a unified 
physical picture of the process of confined quasi-lD scat- 
tering at low energies and can be applied to impurity 
scattering in mesoscopic 2DEG systems and two-body 
collisions. Following the reasoning related to Ea. (|12|) one 
expects to be able to treat non-cylindrical geometries and 
quasi-2D scattering in an analogous way. By computing 
then a few coefficients and functions such as ain and 



Ac (or A„) and only the lowest eigenstates of the confin- 
ing potential, reasonable numerical results are also possi- 
ble. Scattering resonances such as total reflexion, as well 
as weakly localized states can be determined. For the 
particular case of parabolic confinement, the formalism 
presented here can also be used to obtain e.g. confine- 
ment induced two-body weakly bound states. Of partic- 
ular interest is the study of the poles of the scattering 
amplitude which correlate to orbital angular momenta 
beyond s-waves and how non-parabolic geometries can 
affect these unconventional pairings due to couplings to 
the center of mass coordinates. 
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